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Abstract 

Using hydrodynamic approach, it is shown that the properties of a 
marginally stable collisionless stellar disc resemble those of a thermo- 
dynamic system undergoing a gas-liquid phase transition. The max- 
imum in Toomre's stability diagram, which separates gravitationally 
stable and unstable states with respect to axisymmetric perturbations, 
can be treated as a critical point for this transition. Static perturba- 
tions of stellar density are explored and the mean perturbation am- 
plitude is considered as the order parameter of the theory. The disc's 
state is assumed to change as the disc passes through the critical point. 
Since the disc tends to retain hydrostatic equilibrium, structures can 
be formed spontaneously, identifiable with a seed spiral structure. A 
power-law scaling of the order parameter in the vicinity of the critical 
point has been found. The susceptibility and other Landau- Weiss ex- 
ponents similar to those in the Van der Waals theory are calculated. 
The critical behaviour of marginally stable discs at the initial stage 
of their evolution occurs in numerical simulations where snapshots of 
stellar positions reveal stellar splinters and crescents diverging from 
the disc centre. These structures can be a result of the phase tran- 
sition. In numerical simulations, these structures eventually reduce 
to decaying worm-type features because of the 'heating' most likely 
resulting from instability of stellar orbits due to resonances. Under 
favourable conditions the critical behaviour leading to the establish- 
ment of order in a stellar disc can result in the generation of a spiral 
structure. 
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1 Introduction 

Our aim is to clarify whether an ordered state can arise spontaneously in a 
marginally stable stellar disc. We suggest that the marginal state can be a 
critical one with respect to a phase transition. We show that a self-gravitating 
galactic stellar disc is most sensitive to static perturbations of gravitational 
field at a certain length scale when the disc is in a marginal state (where 
chaotic stellar velocity is equal to the local gravitational momentum ttGu/k, 
i.e., the mean gravity force times the epicycle period, where G is Newton's 
gravitational constant, a the surface mass density and k the epicyclic fre- 
quency). Perturbations of stellar density grow with the difference between 
the radial velocity dispersion and its Toomre's characteristic value, and can 
modify the symmetry of the disc when their scale is comparable to the size 
of the stellar epicycle, the natural wavelength of oscillations. 

The susceptibihty of the stellar 'gas' is singular at the top point of the 
well-known Toomre's stability diagram (Toomrc 1964). Moreover, the static 
response to small perturbations of the gravitational potential can be singu- 
lar in the region of instability in the diagram. Therefore internal physical 
equilibrium of a self-gravitating disc requires an analysis deeper than that in 
terms of the standard theory of dynamic instability. 

Guided by physical considerations (Saslaw 1987), we assume that the 
marginal state of the disc is similar to a critical point of the gas-liquid phase 
transition (as a special case of the phase transition of the second kind), and 
the disc's state in unstable region is similar to the condensed state of the van 
der Waals gas. (We stress that this is just a formal similarity; the physics 
of the stellar disc is quite different from that of a gas-liquid system.) Note, 
however, that we consider not the loss of thermodynamic stability as density 
increases, but rather the occurrence of a susceptible state, compressible and 
'elastic' simultaneously. Therefore, the physical properties of perturbations 
in the 'van der Waals region' are analyzed using a nonlinear dependence of 
'pressure' perturbations on stellar density perturbations. We demonstrate 
the possibility of a quasi-static compression of the perturbed stellar 'gas', 
obtain an equation for the order parameter (viz. the amplitude of the static 
perturbations in the stellar spatial density), and investigate the perturbations 
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in the neighbourhood of the critical point in Toomre's diagram. 

In fact, a stellar disc is a collisionless dynamic system. The 'elasticity' 
of such a self-gravitating system results from both the epicyclic motions of 
stars under small perturbations and the centrifugal 'elasticity' of the reference 
(circular) stellar orbits themselves. 

Although critical fluctuations can be important because of an anoma- 
lous behaviour of the susceptibility (Ivannikova, Ivanov & Maksumov 1994, 
1997), we shall restrict ourselves to Landau's approach and consider the pos- 
sibility, for a self-gravitating stellar system, to reach a new local equilibrium 
by developing an ordered state. This approach uses an order parameter as 
a measure of evolution in the disc structure as it undergoes a phase tran- 
sition, but fluctuations are ignored. Below we use the mean amplitude of 
density perturbations in the marginal disc as the order parameter. We do 
not use the method of thermodynamic potentials because of the well-known 
complications in its application to self-gravitating systems (Saslaw 1987). 

As follows Saslaw's arguments (§§29.1, 29.3, 30.3, and 31.1 in Saslaw 
1987), a direct introduction of gravity into the theory of the van der Waals 
gas leads to a poor approximation, especially because of the collective nature 
of the gravitational interaction. 

The method used here combines the mean field approach, similar to the 
mean molecular field theory of Weiss (§9.2 in Balescu 1975), and the ideas 
related to hydrostatic equilibrium in an external field and Landau's theory of 
the gas-liquid phase transition. This differs from the fluctuational approach 
as the perturbations are considered to be ordered in space. It is assumed that 
all the stars are identical, and so each equally interacts with the other stars 
in the disc. This interaction is described by the mean field. This method is 
not free of disadvantages similar to those of the theories of van der Waals, 
Weiss and Landau. 

In Section 2, density perturbation in the pre-critical state is obtained us- 
ing linear collisionless hydrodynamics (theory of density waves, Lin & Shu 
1964) and then its limiting value corresponding to the top point of Toomre's 
diagram is taken. This approach is equivalent to the local gravithermody- 
namics, as shown by Saslaw (1987, §31) for the Jeans criterion as an example. 

But the direct analysis of intrinsic (or proper internal) and forced static 
perturbations (as opposed to fiuctuational ones) at the critical point in terms 
of general thermodynamic relations in Section 3 apparently provides the most 
straightforward explanation of the new phase in the marginally stable disc. 

Our results, albeit schematic, clarify an important question as to how a 
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self-gravitating disc will respond to a static spiral perturbation of gravita- 
tional potential (in the approximation of tight winding) whose wave number 
is equal to the inverse Coriolis radius and how the amplitude of stellar den- 
sity perturbations will vary if the disc is close to the marginal stability (i.e., 
where the mean chaotic stellar velocity is comparable to the mean gravita- 
tional momentum acting on stars that move along the Coriolis circles). We 
assume that the response of the disc is similar to a phase transition of the 
second kind, and we use this assumption to suggest the form of the 'pressure' 
perturbations in the post-critical regime. The marginal state appears to be 
critical with respect to compression driven by the spiral perturbation of the 
gravitational potential, which can be expressed as a growth of the suscepti- 
bility of the stellar density perturbations. We recall that the critical point of 
the gas-liquid phase transition is characterized by the loss of thermodynamic 
stability as density increases. 

The stellar disc in our theory does not have the van der Waals equation 
of state. We found it difficult to obtain any explicit form of the equation of 
state for the stellar disc. However, it is sufficient for our purposes to know 
the 'pressure' perturbation alone. Our main tools are Toomre's stability 
diagram, the density wave theory, and the theory of critical phenomena due 
to van der Waals, Weiss and Landau. Important for the method used are 
Toomre's critical dispersion (as an indicator of criticality) , the mean field 
approach borrowed from the theory of stellar collective phenomena, and the 
ideas related to hydrostatic equilibrium in an external field. 

We further assume that, in the critical region, the disc evolves into a new 
equilibrium state similarly to the van der Waals gas. Therefore we consider 
here static perturbations in the stellar disc. The amplitude of the pertur- 
bations is determined from equilibrium equations in a constant external and 
internal gravitational fields. We recall that Landau considered the thermo- 
dynamic potential for a fixed value of the order parameter determined from 
the minimum condition for the potential. 

The concept of the critical point is modified because of the long-range 
nature of the gravitational force and the absence of stellar collisions (Saslaw 
1987). An important role here is played by Toomre's velocity dispersion, that 
is by the collective local interaction, rather than by pairwise interparticle 
interactions of van der Waals. The 'centrifugal elasticity' also results in 
a nonlinear term in the 'pressure' perturbation. We consider spontaneous 
condensations ('convergence') and rarification ('divergence') of stellar orbits 
in radius (similar to compressible perturbations in the van der Waals gas). 
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but do not include such effects as 'heating'. 



2 The response function, coherence length and 
susceptibihty in the pre-critical region 

In this section we evaluate the response function using linear collisionless 
hydrodynamics including both forced and internal (proper) perturbations in 
the stellar spatial density. 

Consider a nonaxisymmetric perturbation of the stellar surface density 
Cm (the order parameter), whose background value is a, induced by small 
perturbations of gravitational potential ipi = V'i,int + V'l.ext (positive definite), 
where the potential is represented as a sum of internal, self-consistent field 
and external field. In the limit of tight winding, we have (see Appendix A) 

c — ((jj — niQ) 

where Cr is the radial velocity dispersion, Q is the angular velocity of the disc, 
K is the epicyclic frequency, u is the angular frequency of the perturbations, 
and m/r and k* are the azimuthal and total wave numbers, respectively. 
Using V'l.int = '2'T!'G(7m/\k*\, Eq. (P) can be rewritten as 

— [K^~{io-mn)^] = {-clk*^ + 27rGa\k*\)— + ^i,^tk*^. (2) 
a cr 

Assuming that u — mQ = (cf. Saslaw 1987, §31) and dividing by k^, we 
obtain 

Consider the density response (positive definite) and complete the square 
in k* on the left-hand side: 

Clk*^ 2'kGg TT^GV^ TT^GV^ 

1 + — — - H 

ru L-f rv rv rv rv 



1 + 



K 



ttGcx^ 
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where we have taken into account that k* — nGa/c^ = k. 

We can introduce a characteristic length (analogous to the correlation 
length) quantifying the spatial coherence of the perturbations in the mean 
self-consistent field. Two quantities with the dimension of length can be 
identified in the above expression: r^^^.^. = c^/k^ and Ll^^^ = r'l^„cl/{cl — c\). 
In terms of i^corr; Eq. (EI) reduces to 

(Trn _ V'l.ext k*"^ 
^ -^corr + 

This analysis resembles the theory of Lin and Shu of spiral density waves 
in its hydrodynamic form. But there is a crucial distinction. We do not 
address the question of instabihty controlled by the magnitude of the radial 
velocity dispersion. Instead, we consider two types of response in the stellar 
density and its scaling with the velocity dispersion: a spontaneous response 
not affecting the disc energy and a forced response. 

Having accepted Lcorr as the coherence length, we see that the critical 
exponent v is equal to 1/2, as it should in Landau's theory; from the definition 
of z/, Lcorr ~ (c^ — c\)~" . For Cr > ct, 1.6., in the pre-critical range of c^, 
Lcorr is close to the stellar Coriolis radius as obtained by Ivannikova et al. 
(1997) from different arguments. However, the expression obtained here can 
be extrapolated to the critical region. 

Lcorr can also be identified with the spatial extent of the pattern since 
its inverse is the width of the spread in the wave number space, and the 
expression in parentheses on the left-hand side of Eq. (|3]), if equated to zero, 
yields Toomre's dispersion relation for Lin's spiral density waves in the ap- 
proximation of collisionless hydrodynamics with u — niQ = 0. This is the 
equation of the marginal curve shown in Fig. 1. 

Now consider the response in the vicinity of the top point in Toomre's 
diagram. Assuming = in Eq. i.e., k* = k^r = iiGa/c^., we obtain for 
the amplitude of a mode whose scale is close to the critical value k~^: 

(J 
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(4) 



(5) 



Figure 1: Toomre's stability diagram for a thin stellar disc in the hydro- 
dynamic approximation. The region of instability lies below the marginal 
curve. Q is Toomre's parameter, A = 27r/A;, and Acr = At- 

This yields the critical exponent as 7 = 1 for > ct from the definition 
<^m = X''/'i,ext of the susceptibihty x and the definition of 7, x ~ ~ 
c^)'"' . As expected, the correlation length is related to the susceptibility (Ma 
1976, Patashinsky & Pokrovsky 1975). The spontaneous value of the order 
parameter cTm/c" in the post-critical state is also related to L^on (Patashinsky 
& Pokrovsky 1975), i.e., is determined not only by the direct local correlation, 
but also by the global correlation chain whose scale grows as the critical point 
Cr = Ct is approached. 

There are further implications of Eq. For Cr < ct, this equation is 
inapplicable if k tends to zero remaining on the marginal curve because then 
Fredholm's alternative for the inhomogeneous problem is violated. Physi- 
cally this means that matter motions in the disc must be allowed, for which 
nonlinear analysis would be required. 

Two types of the critical behaviour of the disc can occur. Firstly, the disc 
can become dynamically unstable, leading to processes similar to those stud- 
ied by theory of synergetics. Secondly, a phase transition can occur, similar 
to the gas-liquid phase transition at a critical point. Here we investigate 
the latter posibility. 
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3 The order parameter, susceptibility and crit- 
ical exponents in the critical and post-critical 
states 

In the previous section and Appendix, the disc was treated in the hnearized 
WKB approximation as an isothermal disc. In this section a different, non- 
linear expression for the 'pressure' perturbation must be used. The reason 
is that the properties of the disc matter change at = ct. The disc has 
different properties in the pre-critical state (c^ > ct), in the marginal one 
{cr = Ct) and in the post-critical state (cr < ct, \cr — ct| << ct). In 
the latter state, the mean field V^i^mt is no longer destroyed by the peculiar 
motions of the stars and, presumably, affects the disc properties. One can 
use the linear response function when > ct as in Section 2. However, 
the linear approximation (linear perturbation analysis) must be abandoned 
when we consider in this section the general equilibrium state for the marginal 
and post-critical regimes. We consider perturbations of a small, but finite 
amplitude. 

Let us return to Eq. (|2]) . As follows from the linear theory of density waves 
and from the above analysis, the most unstable (critical), as well as the most 
susceptible wave number is given by fc^uns — in the marginal state 

with u — mQ = and Q = 1, where fcy = K^/(27rG(j), Q = Cf/cx, and 
Ct = (vrGcr)//? (see Fig. 1). For this mode, Eq. ([2]) takes the form 

a a 

Using fcm.uns = i^/ct, canceling and after some algebra, Eq. (EJ) can be 
reduced to ^ 

— {cl - Ct) = ^l.ext • (7) 

a 

Equation ([7]) correctly describes the susceptibility of the critical mode 
for Cr > Ct but does not yield spontaneous phase transition (at V^i^ext = 0). 
Furthermore, it does not describe the susceptibility in the post-critical state 
Cr < Ct since it has been obtained without taking into account changes in 
the disc's state. 

Now we leave the mode analysis for the general equilibrium relations. We 
recall that in fact a stellar disc is a coUisionless dynamic system. 

The relation between 'pressure' perturbation Ap and gravitational po- 
tential function ipi applicable to static perturbations over half the Toomre 
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wavelength (i.e., the critical scale) has the form 

^-*,=0. (8) 

a 

Such perturbations can be called frozen. They satisfy the general hydrostatic 
equilibrium equation for a system in an external field known from statistical 
physics (Landau & Lifshitz 1976, §25): 

fio{P,T) - ^(r) = const , 

where /zq and \1/ are the specific chemical potential (without the field) and 
field potential, respectively, P is the 'pressure', and T is the temperature of 
the system. \1/ is reduced by the centrifugal potential (cf. Eq. (pHj) of the 
Appendix A that describes hydrostatic equilibrium of the disc in the pre- 
critical state). 

The 'pressure' perturbation near the critical point (i.e., the top of 
Toomre's diagram) can be written in a form similar to that used by Landau 
in his theory of van der Waals' critical point (Landau & Lifshitz 1976, §152): 

c;^ + 44^. (9) 



.2 

a a a 



The applicability of this expression to stellar discs is our assumption, justified 
a posteriori by the results obtained. This expression follows from general 
arguments presented, e.g., by Landau & Lifshitz (1976, §§144, 152). For 
static perturbations, this expression results in one real and two imaginary 
solutions for am when > ct and three real solutions when c^. < ct- Thus, 
the expression (jH]) leads to a reasonable physical behaviour similar to that of 
a gas-liquid system, in the sense that the stellar disc maintains hydrostatic 
equilibrium near its marginal state. The expression ([9]) simply implies the 
expansion of Ap in power of am, supposing the inflection (bending) point 
for equation of state function of a stellar disk p{a) (in a marginal state), as 
Ap = {dp/da)am + {^/n\){d^p/da^)a'^. That is, the disc's state is assumed 
to change as the disc passes through the critical point, as we pointed out in 
Abstract, the inflection (bending) point is appeared . 

The term nonlinear in am in the expression can be understood in 
the followimg semi-quantitative way. A density perturbation am leads to 
an increase in stellar velocities on their orbits by Av due to an increase 
in the gravitational force. For static perturbations, the velocity incere- 
ment can be found from centrifugal balance as QAv = irGam- Thus, Av = 
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-nGoml^ ~ I'nGo^j K. The resulting 'pressure' increment is given by 
(Af)^o"m ~ 4(7rGcr/K)^(crm/cr)^o"m, which is just the nonhnear term of the 
expression ([9]). This term represents an additional 'elasticity' (apart from 
that of the epicyclic motion, represented by c^) that arises from the stability 
of the circular orbits, associated with the centrifugal force, because the orbits 
are assumed to sustain gravitational perturbations for c^. < ct- 
We note that 



^l.int 

which results, using Eq. 

a 



^l,mt + ^1, 



ext 5 



m' corr 



Ct 



K 



a 



in the following equation for the order parameter 



- Ct + 4ct 



a 



^1. 



ext 



(10) 



Equation ffTUl) has a standard form for an equation describing phase equi- 
librium in Landau's theory of the phase transition of the second kind (Landau 
& Lifshitz 1976). This equation yields the ratio controlling the susceptibility 
d<Jm/d4'i,ext in the pre-critical state {cr > ct), 



dar, 



X 



ext 



^2 ' 



(11) 



together with expressions for the dependence of the order parameter on c^—c^ 
(if evaluated at V'l.ext = 0, see Eq. f lT^ below) and on the external field (if 
evaluated at c^ = ct; note that the result remains non- singular ) , 
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Thus, there are three critical exponents, 7 = 1, (3 = 1/2 and 6 = 3. These 
are not independent but related by Widom's ratio (Landau & Lifshitz 1976, 

§148),/35 = /3 + 7. 

For Cr < Ct, the susceptibility has to be found for a non- zero value of 
am- For this purpose the bracket in Eq. f lTOj) should be put equal to zero. 
For Cr < Ct we obtain a real solution of the type am/cr ~ (ct — c^)^ which is 
given by 



a 



cl 



(12) 
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For Cr > Ct, Eq. fll2l) formally gives purely imaginary density perturba- 
tions, implying that static density perturbations of finite amplitude cannot 
exist in this case. 

Differentiating Eq. ffTOl) and using Eq. f lT2|) . we find the susceptibility in 
the post-critical case as 



dam _ _ 0^ 

a 



m - 4) 



(13) 



Comparing Eqs ffTTj) and f lT3|) . we see that the susceptibility suffers a jump 
at Cr = Ct as it should at a phase transition. 

Finally, the exponents follow as/3 = l/2, 7 = 1, 6 = 3, and u = 1/2. The 
scalings of Eqs ( 1TT|) -(1T31) agree with the results of Landau- Weiss's standard 
theory (Landau & Lifshitz 1976). 



4 Conclusions 

In the spirit of Landau's theory of the van der Waals critical point and the 
generic nature of critical phenomena, we have shown the possibility of critical 
behaviour of galactic stellar discs. The disc undergoes a phase transition 
when the stellar velocity dispersion is equal to Toomre's critical value, 
Ct. It appears that the transition leads to the formation of a fundamental 
spiral structure. We have shown that well-known features of such phase 
transitions result in an essentially new mechanism for the spiral structure 
formation. Similarly to a thermodynamic system undergoing the gas-liquid 
phase transition at a critical point, a marginally stable coUisionless stellar 
disc can, while retaining (or tending to retain) its hydrostatic equilibrium, 
exhibit a critical behaviour that leads to its primary spatial inhomogeneity 
which can further develop into a spiral structure. This is the main point 
of our approach which takes advantage of the theory of critical phenomena; 
linear results are a part of our analysis, but they are not sufficient for it. 

Numerical experiments of Hohl (1975, 1975a) demonstrate that a stellar 
disc embedded in a bulky halo appears to be rather changeable ai Q = Q^- 
Processes different from the bar instability 'thermalize' the disc quickly (prac- 
tically instantaneously), heating it up to Q ~ (2-3)Qcr- At the initial stage 
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of the process, the disc is populated by a set of stellar splinters and cres- 
cents, and a few large-scale spiral branches develop later. These structures 
eventually reduce to decaying worm-type features because of the 'heating' 
most likely resulting from the instability of individual stellar orbits due to 
resonances. This demonstrates that both dynamical instabilities and fluctu- 
ations (critical and, for example, induced by embedding bulky masses) can 
be superimposed on the process wherein an order emerges. This scenario 
is also confirmed by numerical calculations of critical dynamics in a model 
one-dimensional self-gravitating system in a marginal state (Ivanov 1992, 
1993). 

'Hot' discs cannot have quasistatic structures, and this fact corroborates 
our theory of phase transitions. Our theory is only applicable to early stages 
of the simulations discussed above before instabilities have developed. There 
are many effects that can halt the ordering process, so it can occur only in 
some galaxies. 

We have considered only compressional perturbations. But the last term 
in the expression for Ap/a can be compared with a change in the stellar 
velocity driven by gravitational perturbations during the phase transition, 
namely with the drift correction to the chaotic stellar velocity, whose az- 
imuthal component is given by \K,~^dipi^i^t/dr\ and the radial component by 
\{K,r)^^dil!i^int/df\ (Ivannikova & Maksumov 1994). The square of the drift 
correction to the velocity is given by |/t~"^9'?/'i,int/f^^P + |(^'^)~^c^'^i,int/f^V^P- 
Using T/'i^int = 27iGam/\k*\, k^^^^^^ = '^/ct, \i^~^dilJi^int/dr\ = 2cT(cr„/a) cos (5 
and {Kr)~^dipi^int/dip\ = 2cT{(Jm/o') sin 5, the drift correction squared can be 
expressed as 4cT(crm/o")^. This estimates confirms the possibility of transla- 
tional and shearing effects in addition to the compressional one. 
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5 Appendix A: The derivation of Eq. (1) for 
density perturbations in a stellar disc. The 
drift perturbations 

Linearized hydrodynamic equations for the density perturbation am and the 
radial and azimuthal bulk velocities u and v are given by (Rohlfs 1977, §5.5) 



darn ^ darn 1 

ot Oif r 



-{ran) + -{av) 







du ^ du 
at dip 

dv dv 



dam ^ dipi 
a dr dr 



1 / c; dar 



(14) 

(15) 
(16) 



r \ a dip dip J 

Assuming that the perturbed quantities (including ipi) are proportional to 
exp[i{ujt — rriip + kr)] and reserving for amplitudes the notation used for the 
respective perturbed quantities, we transform these equations first similarly 
to Lin & Shu (1964) who considered tightly wound perturbations. Equations 
(111, ([151) and ([161) reduce to 



(17) 

(18) 
(19) 



i(uj — mil)am + iaku + —u — iak^v = 

dr 

kc^ 

i — -am + i{<^ — mQ)u — 2flv = ikipi , 
a 

i{io — m\l)v H — —u = I -ar 

^ ' 2VL a 



ik^il) 



Equation f[T9|) yields for v: 



K 



k^c^ am/o- - k^ipi 



2il(u — mil) ' u — mVt 
Using this expression in Eq. (ITS|) . we find u and v as 



u 



oj — niQ 



X 



-kicj^ 
a 



^1 
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{uj — mVtY — ft;^ 
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' cr 



2Vt 



k,. 



a 



^1 



Using cl 



UJ — mVt 

(KV2fi)c2, we find from Eq. ([17]): 



(J 



1 



1 dhiVt a„ 



2 dlnr a 



+ zizfo; — miZ) — - — 
dr 

■ 11, dfl . 

+ ikk^ioo — mil) r—— ipi 
dr 



K 

2n 



2 '^m 



^1 



—i{uj — mn)' 



din a 

dr 



k ( d— - tpi 
a 



Neglecting terms with spatial derivatives of a and Q, we obtain Eq. ([T]). 

Now consider slow, drift modes. They are basically different from the 
tightly wound ones and can only be nonaxisymmetric. Their susceptibility 
is easier to consider in the limit d/dt + Qd/dip = 0. The system flT^ - ffTUj) 
then takes the form 

(rcru) + —((Tv) =0 , (20) 



dr 

2nv 



K 



a dr 

1 



c; <9cr„ 



dr ' 



(21) 

a dip dip ) ' ^^^"^ 

In the limit m, f — t- we obtain an expression used to determine the suscep- 
tibility of the drift modes: 



2 '-^m 



a 



-^1 
14 



. 



(23) 



Using ijji = V'l.int + V'l.ext and V'l.int = 2TcGam/\k*\ with k* = Ak^, k* = 
{k'^ + k'^) , we obtain Eq. (|TT|) for the susceptibihty x- 
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